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I will review the measurements of the 6-quark mass performed by ALEPH and DELPHI. A large set of ob- 
servables has been used together with detailed studies on jet algorithms. Very clear effects due to the 6-quark 
mass running are observed. Comparing with the determinations at the T(4S) threshold, the measurements of the 
running 6-quark mass at the Z pole are consistent with the predicted evolution from QCD. 



1. Introduction 

The 6-quark mass is one of the fundamental pa- 
rameters of the QCD Lagrangian. However, due 
to confinement, quarks do not appear as free par- 
ticles and therefore the definition of their mass is 
ambiguous. In fact, quark masses can either be 
defined as for free particles (such as the leptons) 
as the position of the pole of the propagator, or 
they can be interpreted as effective coupling con- 
stants in the Lagrangian. In the former definition 
the mass is called "pole mass" and does not run 
with energy; in the latter the mass is called "run- 
ning mass" ; as such it is a scale and renormaliza- 
tion scheme dependent quantity. 

The 6-quark mass is important for phenomenol- 
ogy since it enters many theoretical predictions 
of physical quantities such as CKM matrix el- 
ements, &-hadron inclusive semilcptonic decays, 
total widths. 

Most of the 6-quark mass determinations have 
been extracted at the T(4S') threshold, from the 
T bound states with QCD sum rules [|l| and QCD 
lattice calculations 0. It is therefore interesting 
to measure this parameter at higher scales such 
as the one offered by LEP. 

At LEP, the possibility to test the running of 
the fe-quark mass within the framework of the 
perturbative QCD had not been considered until 
very recently. The reason is that the effects of the 
mass rapidly become very small with increasing 
energy for many observables, since they are pro- 



portional to to^/M| (O (0.1%)), for instance in 
the total cross section. For other quantities, such 
as jet rates, the effects are of the order of few per- 
cent due to the suppression of the radiation off b 
quarks. The advent of theoretical calculations at 
next-to-leading order (NLO) in QCD perturba- 
tion theory allowed for the determination of the 
&-quark mass [|| . 

In this article I will review the measurements 
done by ALEPH 0] using a large set of event 
shape observables and by DELPHI P,p| using the 
three-jet rate. The different sensitivity to the 
hadronization corrections and next-to-leading or- 
der corrections can be used to investigate the in- 
terpretation of the results. 

2. Analysis Method 

The method to extract the 6-quark mass is 
based on measuring the ratio of any infrared safe 
observable O computed for b and uds (referred as 
q) induced events and assuming a^ universality. 
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The following variables have been studied by 
ALEPH: 

• The rate of three-jet events, where jets are 
reconstructed using the DURHAM cluster- 
ing algorithm with i/cut = 0.02. 

• The first and the second moments of 
the event shape variables Thrust T, C- 
parameter C, j/3. Total and Wide Jet Broad- 
ening, Bt and Bw- The definitions of the 



observables can be found, e.g., in |H,g| and 
references therein. 

DELPHI has used the three and the four jet 
event rates. 

The measured ratio R'^J^'^^ of the observables in 
^-tagged events and in uds-tagged events can be 
related to the quantities at the parton level via 
the following formula 
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{ObHbDbTbbPbb + OcH,D,TbcPbc + 
OgHgDgTbqVbg) / {Ob Hb D bTgbV gb+ 

OcH,D,Tq,Vqc+ OqHqDlTqqVqq) .(2) 



Here Hx,Dx^Tyx are the corrections due to 
hadronization, detector effects and tagging, re- 
spectively, and Vyx are the purities of the tagged 
sample, where x is the true flavour and y is the 
tagged one. The tagging corrections take into ac- 
count biases introduced by the flavour tag. In 
case of ALEPH the ratio is computed with re- 
spect to an unbiased sample of all flavours, which 
means that Vqx is equal to the ratio of the partial 
width of the Z to X quarks and the total hadronic 
width, and Tq^ — 1. 

All the correction factors and purities are ob- 
tained from Monte Carlo (MC) simulations. Be- 
cause mainly ratios of corrections are involved, 
some systematic uncertainties cancel. 

^T i^ extracted from the relationship (|^) and 
finally corrected for the contribution of anoma- 
lous triangle diagrams jl^l, in order to relate 
B^'l' to R^'^ for which the perturbative calcula- 
tions have been performed. 

The general form of the NLO prediction for 
■^ll ^^ ^ function of the runninng 6-quark mass 
nib is of the form 
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The ratio can also be expressed in terms of the 
pole mass Mb 
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The two predictions are equivalent at this or- 
der. The coefficient functions 6o,i for the two 



schemes can be related to each other by express- 
ing the pole mass in terms of the running mass 
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The coefficient functions 6o and 6i have been 
computed for the ratio of three and four jet rates 
p,[ll|. The three jet rate ratio allows a NLO pre- 
diction, while for the four jet rate only a LO pre- 
diction is available. For the other variables used 
by ALEPH the predictions have been obtained us- 
ing the MC generators ZBB4|ll and EVENT |l|] 
which are correct to NLO. The hadronization cor- 
rections have been evaluated for all variables by 
computing the relevant observables at parton and 
at hadron level. 

3. ALEPH analysis 

In the ALEPH analysis data taken at the peak 
of the Z resonance from 1991 to 1995 are used. 
A standard hadronic event selection [g[ is ap- 
plied, which is based on charged particles. A cut 
|cos6't| < 0.7 is imposed, where 9t is the po- 
lar angle of the thrust axis, computed from all 
charged and neutral particles as obtained from 
the energy-flow algorithm |I4| . This requirement 
ensures that the events are well contained within 
the vertex detector (VDET) acceptance. After 
the selection, a sample of 2.3 million hadronic 
events remains for further analysis, with about 
0.3% of non-hadronic background, mainly com- 
ing from T~^T~ events. 

An important ingredient of the analysis is the 
choice of the Monte Carlo simulation, which is 
based on on JETSET 7 . 4 parton shower model 
plus string fragmentation fls]. The production 
rates, decay modes and lifetimes of heavy hadrons 
are adjusted to agree with recent measurements, 
while heavy quarks are fragmented using the Pe- 
terson et al. model |lq| . 

The observables described in Section g are com- 
puted using only charged tracks. 

In the analysis only b events are tagged. The 
selection is based on the long b hadron lifetime 
thanks to the great precision of the VDET, and 
follows closely the approach described in detail 
in Ref. [|^. The actual selection of b events is 



Table 1 

Leading (LO) and next-to-leading (NLO) order contributions to 1 — R^^'^ , for the running mass (run) 
and the pole mass (pol) schemes. The contributions are evaluated for a running (pole) mass of 3 (5) 
GeV/c^. The strong coupling as{Mz) is set to 0.119. The values are given for the three-jet rate (DURHAM 
scheme) and for the event-shape variables thrust T, C parameter, the transition value 2/3 for three to two 
jets (DURHAM scheme), and the total and wide jet broadenings {Bt and Bw)- The indices indicate the 
first or second moment of the event shape variable. 
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0.056 


-0.008 


Ti 


0.036 
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-0.007 


T2 


0.017 


0.032 


0.043 


0.036 


Cx 


0.044 


0.022 


0.091 


-0.011 


C2 


0.021 


0.039 


0.052 


0.043 


J/3i 


0.032 


0.007 


0.071 


-0.021 


2/32 


0.015 


0.003 


0.032 


-0.007 


Bti 


0.117 


0.006 


0.188 


-0.074 


Bt2 


0.036 


0.112 


0.080 


0.123 


Bwi 


0.117 


-0.085 


0.188 


-0.183 


Bw2 


0.036 


0.016 


0.080 


-0.013 



obtained from a cut on the distribution of the 
confidence level Puds that all tracks of the event 
come from the main vertex. Very good agreement 
between data and MC is observed. On b events 
the selection has an efHciency of 80.5% and a pu- 
rity oiVbb = 83.1%. The backgrounds amount to 
Vbc = 13.5% and Vbg = 3.4%. 

3.1. Choice of variables 

A list of the leading (LO) and next-to-leading 
order (NLO) contributions to 1 — ii^d'^* is given in 
Table |l] for all variables, both for the running and 
the pole mass schemes. They have been evaluated 
for a 6-quark mass of 3 GeV/c^ in the former and 
5 GeV/c^ in the latter case. It is found that for 
some observables such as the jet broadening vari- 
ables the mass effect is rather large. However, the 
NLO corrections can also be sizeable, as in the 
case of the second moment of thrust and the first 
moment of the wide jet broadening. For such vari- 
ables it would definitely be necessary to also com- 
pute the NNLO contributions in order to obtain a 



reliable perturbative prediction. Because of these 
observations only those variables are considered 
for which in both schemes the NLO contribution 
is clearly smaller than the LO term. 

3.2. Hadronization corrections 

The perturbative predictions are corrected for 
hadronization effects by computing the relevant 
observables at parton and at hadron level, in- 
cluding final state photon radiation off quarks. 
Several Monte Carlo models based on the parton 
shower approach plus subsequent string or cluster 
fragmentation are employed for this purpose. For 
the nominal analysis, the same generator as for 
the full simulation is used. 

The ratios of hadronization corrections, -ffb/g, 
are listed in Table ||. The corrections are rather 
sizeable for almost all the observables; in most 
cases they are of the same size as or larger than 
the expected mass effect. Only the three-jet rate 
and the first two moments of the 2/3 distribution 
have corrections at the 2% level or below. For 



all the other event-shape variables the corrections 
are of the order of 10% or even larger. It has been 
found that the deviation from unity is almost en- 
tirely due to b hadron decays, which change the 
distributions mainly in the two-jet region. As 
can be observed from Table g, the same correc- 
tions, computed taking only hadrons which stem 
directly from the string before any decay, come 
close to unity within one or two percent. This is 
in agreement with expectations from recent calcu- 
lations of nonperturbative power-law corrections 
to moments of event-shape distributions for mas- 
sive quarks [ p^ . 

By requiring that the hadronization corrections 
do not exceed 2% two variables remain: the first 
moment of 2/3 and the three jet ratio. 



3.3. Measurement of R^^'^ 
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The results for the pcrturbative ratio i?j,^ 
listed in Table || after having applied the correc- 
tion procedure defined in equation [^ and correct- 
ing for anomalous triangle diagrams. 

The dependence on the resolution parameter 
ycut of the pcrturbative ratio R^^'^ for the three- 
jet rate is indicated in Fig. |[ 

The data are compared to the predictions of 
the parton shower models of PYTHIA 6. 1, which 
is based on JETSET, and HERWIG 6.1 [|l|]. The 
MC models are in reasonable agreement with the 
measurement at large j/cut values, where the sta- 
tistical uncertainty is large, however. At interme- 
diate resolution parameters these models predict 
a lower ratio than observed. 

In addition, in Fig. n the next-to-leading or- 
der pcrturbative QCD predictions for two differ- 
ent values of the running 6-quark mass in the MS 
scheme are given. The data clearly favour a mass 
close to 3 GeV/c^ 

3.4. Systematic errors 

Systematic uncertainties on the mass extrac- 
tion can be divided into three categories: the 
ones coming from uncertainties on the tagging bi- 
ases, from the hadronization and from the missing 
higher order corrections. 

The first ones have been evaluated by varying 
the gluon splitting rate which influences the pu- 
rity of the b tag and by studying the variation of 
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Figure 1 . Comparison of the j/cut dependence of 
the measured ratio R^'^ for the three-jet rate to 
the predictions of parton shower models as well as 
next-to-leading order (NLO) perturbative QCD 
for two different values of the 6-quark mass in the 
MS scheme. The errors are statistical only. 



the results due to the inadequacy of the simula- 
tion of the tracking detectors. 

The uncertainties from the modelling of the 
hadronization are typically evaluated by comput- 
ing the hadronization corrections with different 
MC generators. As has been shown in Section 
3.2, the b hadron decays have a large impact on 
the size of the hadronization corrections. Because 
in the tuning of the MC hadronic final states are 
analyzed after all decays, differences in the de- 
scription of hadron decays can lead to differences 
in the tuned fragmentation parameters. In order 
to assess an uncertainty related to the string frag- 
mentation parameters and subsequent decays, the 



Table 2 

ALEPH results for i?^" with statistical errors, including the uncertainty from the MC statistics. The 
results are obtained with the hadronization corrections -ffh/g as predicted by the MC. Also given are the 
corrections HJ^)"^, computed from hadrons directly originating from the string, before any decays {nod=no 
decays) . 






RT 
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Rs 


0.974 ±0.005 


0.980 


0.987 


Ti 


0.910 ±0.002 


1.134 


1.003 


Ci 


0.894 ±0.002 


1.169 


1.009 


ysi 


0.955 ±0.005 


1.023 


0.987 


J/32 


0.979 ±0.010 


0.984 


0.984 


Bt, 


0.831 ±0.001 


1.306 


1.025 


Bw2 


0.929 ±0.003 


1.088 


0.986 



variation in the final result when using two differ- 
ent b hadron decay tables for the hadronization 
corrections is taken as a systematic uncertainty. 

Another quantity relevant for this analysis is 
the b fragmentation function, which describes 
the energy fraction transferred to the b hadrons 
during the fragmentation process. Three differ- 
ent sets of different fragmentation functions were 
used to assess the systematic uncertainty due to 
this source. Finally, in order to study purely the 
difference between string and cluster fragmenta- 
tion, the variation in the hadronization correc- 
tions for light quarks (uds) only has been propa- 
gated into an uncertainty on the final result. 

Missing higher order corrections are estimated 
by extracting firstly the pole mass from the per- 
turbative expression in the pole mass scheme, 
then translating that result into a running quark 
mass mii{mi,) at the fe-mass scale, and finally run- 
ning this mass up to the Mz scale using the renor- 
malization group equation. In addition, the ef- 
fects of uncalculated higher order terms can be es- 
timated by a change in the renormalization scale. 
For the central value /i = Mz is employed, and 
the systematic error is taken to be half of the 
range of mass values (average of running and pole 
mass schemes) found when varying /i from O.lMz 
to 2Afz. The behaviour of the extracted mass in 
the two schemes as a function of the renormal- 
ization scale /i is illustrated in Fig. H for the first 



moment of the 2/3 variable. 

3.5. Results for mb(Mz) 

Based on the predictions obtained above, the 
running 6-quark mass is determined from the 
measured ratio i?|^" for two observables, the first 
moment of the j/3 distribution and the three-jet 
rate. 

As stated in Section 0, the functional depen- 
dence of the coefficients 60,1 for the moments of 
the event shape variables are computed by the 
ratio of the cross sections for b and d quarks us- 
ing the MC generators ZBB4 [Q and EVENT 
PJ] . ZBB4 allows for the integration of the fully 
differential NLO matrix elements including mass 
effects, whereas EVENT contains the massless ex- 
pressions 0|. 

In Fig. Wan example of these fits is shown for 
the first moment of the 2/3 distribution. With 
these parametrisations it is possible to estimate 
the actual quark mass effects in leading order and 
its next-to-leading order corrections. In Table || 
a list of these LO and NLO contributions is given 
for all variables for both the running and the pole 
mass schemes. 

The results are listed in Table |^ together 
with the statistical and systematic uncertainties, 
which have been propagated from the correspond- 
ing uncertainties on R^^'^ . The result from the 
first moment of the 2/3 observable is quoted as 
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Figure 2. Renormalization scale dependence for 
the extracted 6-quark mass in the running and 
pole mass scheme, for the observable y^-^. The 
error bar spans the range of the scheme uncer- 
tainty. The full line indicates the renormalization 
scale dependence of the average computed from 
the running and pole mass schemes. 



Figure 3. Parametrisations of the 6-quark mass 
dependence for the LO (top) and NLO (bot- 
tom) contributions to R^^'^ in the running mass 
scheme for the first moment of the j/3 distribu- 
tion. The points indicate the result of the MC 
integration, and the full line the parametrisa- 
tions. The NLO contributions are evaluated using 
as(Mz)- 0.119. 



the final 6-quark mass value, due to its smaller 
hadronization and systematics errors. The the- 
oretical uncertainty is estimated as described in 
Section 3^ by evaluating the impact of the un- 
certainty on the strong coupling constant and the 
renormalization scale variation. Because the pole 
and the running mass scheme are equivalent at 
NLO, the average of the values found for the two 
schemes is quoted as the final result, and half of 
the difference is taken as an additional theoretical 
systematic uncertainty due to the scheme ambi- 
guity. This leads to a measurement of the 6-quark 



mass from ALEPH of 

mb{Mz) ^ [3.27 ± 0.22(stat) ± 0.22(exp) 
±0.38(had) ± 0.16(theo)] GeV/c^. 

4. DELPHI analysis 



(6) 



The first measurement of the running 5-quark 
mass at scales of the Z mass has been performed 
by DELPHI, using 1992 to 1994 data and the 
three jet rate ratio using the Durham algorithm 
for jet clustering Q. Effects of the running of the 
6-quark mass have also been measured in the four 



Table 3 

ALEPH measured &-quark mass mb{Mz) in the MS scheme with statistical (stat), experimental (exp) 
and hadronization (had) uncertainties. In the last column the measured pole mass Mb is listed. 

"O mb{Mz) [GeV/c^] ±(stat) ±(exp) ±(had) Ah [GeV/c^] 
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2/3i 



2.76 
3.34 



0.28 
0.22 



0.28 
0.22 



0.62 
0.38 



3.65 
4.73 



jet rate ratio as reported in |20|. 

A new analysis has been presented at this con- 
ference which is based on the 1994 and 1995 
data sample, which amounts to about 1.3 Million 
hadronic Z decays. A better b tagging; algorithm 
was developed, which is described in [g2[. The ef- 
ficiency to tag b events is about 53% with a 85% 
purity. 

The Monte Carlo model used is based on 
PYTHIA, with the DELPHI tuning. 

The study was performed using the DURHAM and 
CAMBRIDGE^ jet algorithms. The latter is found 
to minimize the NLO corrections to the three jet 
rate for the running 6-quark mass [El| . 

4.1. Measurement of R^^'^ 

The three jet ratio is shown in Figure^ for both 
jet algorithms. As can be seen the LO calculation 
did not provide good description of the data and 
NLO corrections are indeed needed. In the case 
of DURHAM the data are closer to the running b- 
quark mass calculations, despite the NLO correc- 
tions are larger with respect to the pole mass. For 
the CAMBRIDGE algorithm data stiU prefer the run- 
ning 5-quark mass calculations, but here the NLO 
corrections are in fact smaller for this scheme. 

4.2. Systematic errors 

Hadronization errors were derived from com- 
parisons with different Monte Carlo models tuned 
to describe the data. Two different sources of er- 
rors were considered: amod which describes the 
difference of the hadronization corrections be- 
tween PYTHIA and HERWIG, and atun which ac- 
counts for the possible variation of the main frag- 
mentation parameters in PYTHIA. The first uncer- 
tainty takes into account both the different frag- 
mentation scheme (string or cluster) and the dif- 



ferent b hadron decay tables. In the case of crtun 
the error was evaluated by changing the most 
relevant fragmentation parameters of the string 
model in PYTHIA by ±2(t from their central tuned 
values and assuming they are uncorrelated. The 
two above errors were combined in quadrature to 
estimate the hadronization uncertainty. The in- 
terval of validity for the jet resolution parameter 
Vcut to perform the measurement was chosen to 
have small hadronization error (below 5%) and 
smooth Tjcut dependence. Figure |5| shows the to- 
tal hadronization correction uncertainty for the 
three jet rate for the two jet algorithms. 

The uncertainty coming from the tagging is 
taken into account by varying the purities by 1%. 

The errors coming from missing higher or- 
ders are estimated by varying the renormalization 
scale fj, = Mz from 0.5 to 2. The ambiguity due 
to the definition of the three jet rate in terms of 
the running or the pole fe-quark mass has been 



evaluated as explained in Section 3.4 



4.3. Extraction of the 5-quark mass 

The values for the running 6-quark mass ob- 
tained with the DURHAM and CAMBRIDGE algo- 
rithms are shown in Table ^, for i/cut = 0.02 and 
0.005, respectively. 

The three jet rate measured with the 
CAMBRIDGE algorithm is in agreement with both 
LO and NLO calculations in terms of the run- 
ning 6-quark mass. For this reason the mass is 
extracted from this algorithm using the NLO pre- 
diction. The value of the running 6-quark mass 
from DELPHI is therefore 



mb(Mz) = [2.61 ± 0.18(stat)l°:l^(frag.) 
±0.18(tag.) ± 0.07(theo.)] GeV/c^. 
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Figure 4. Corrected data values of the perturbative three jet ratio using DURHAM (left) and CAMBRIDGE 
(right) jet algorithms for DELPHI, compared with the theoretical predictions from B at LO and NLO 
in terms of the pole mass Mb = 4.6 GeV/c^ (dashed lines) and in terms of the running mass mb{M'z) = 
2.8 GeV/c^ (solid lines). 



5. Conclusions 

ALEPH has studied a large set of observables 
to extract the 6-quark mass at the Z peak. The 
first moment of ys has the smallest hadronization 
and systematic errors. DELPHI used instead the 
three jet ratio using the CAMBRIDGE jet algorithm 
which is found to minimize the NLO corrections 
for the running fe-quark mass. 

Clear effects due to the running of the 6-quark 
mass are visible and confirm the expectation from 
QCD. The results from ALEPH, DELPHI and 
a measurement from Brandeburg et al. [g^ are 
shown in Fig. g, together with the average of the 
low energy measurements for mt,{mb) p4| , which 
is evolved up to the Z mass scale using a two-loop 
evolution equation with as= 0.119 ± 0.003. 

The main error on the measurement comes 
from the uncertainties affecting the hadronization 
corrections. A big fraction of it is coming from 



b hadron decays as has been estimated in great 
detail by ALEPH. 

Despite of the fact that the 6-quark mass mea- 
sured from the three jet ratio from ALEPH and 
DELPHI shows a good agreement. 

The mass extracted from the first moment of 
2/3 seems to be larger and barely compatible with 
the jet ratio measurements. 

The spread in the central value and the actual 
size of the errors indicate that the limiting preci- 
sion of the measurement of the 6-quark mass at 
the Z mass scale is of the order of 500 MeV/c^. 
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Table 4 

DELPHI measured 6-quark mass mb{AIz) in the MS scheme with statistical (stat), hadronization (frag.), 
b tagging (b-tag) and theoretical (theo) uncertainties. The first line contains the results of the previous 
published analysis. 
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Algorithm mb(Mz) [GcV/c^] ±(stat) ±(frag) ± (b-tag) ±(theo) Data 
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Figure 5. Evolution of the total hadronization er- 
ror on the three jet ratio for the DELPHI analysis 
with the resolution parameter j/c- The error cor- 
responds to two times the sum in quadrature of 
atun and a modi as described in the text. 
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Figure 6. Comparison of the ALEPH and DELPHI result for TOf,(A/z) with the world average value of 
low-energy measurements for ■mi,{rni,) p4[|, which is evolved up to the Mz scale using a two- loop evolution 
equation with as(Afz) = 0.119 ±0.003. Also shown is the measurements by Brandenburg et al. pJ. The 
inner error bars indicate the quadratic sum of the statistical and experimental uncertainties. The three 
points at the Z pole are separated for clarity. 
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